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DETECTING  LINE-ROAD  AND  ROAD-INTERSECTION 
PATTERNS  AT  VARIOUS  ANGLES 


INTRODUCTION 

Since  1979,  the  U.S.  Army  Engineer  Topographic  Laboratories  (ETL)  has 
been  investigating  discrete  function  technology,  particularly  Walsh  functions,  to  ex¬ 
tract  and  recognize  a  selected  set  of  manmade  cartographic  features,  such  as  road  inter¬ 
sections,  straight-line  roads,  and  rectangular  objects  from  aerial  photographs.1  2  3  4 

The  Walsh  transform  was  chosen  because  of  its  simplicity  (Walsh  functions 
are  only  two-valued,  +1  or  -1),  resulting  in  simple  implementation,  and  because 
Walsh  functions  conform  to  the  selected  set  of  cartographic  features. 

The  technique  was  investigated  by  (1)  using  a  32-  by  32-element,  solid 
state  sensor  array  to  convert  imagery  into  an  electronic  signal,  which  was  processed 
in  a  minicomputer  to  yield  the  Walsh  transform  of  the  images;5  6 '7  and  (2)  using 
a  prototype  image  spectrum  analyzer  (PISA)  with  a  large  size  plasma  discharge 
device  (8.5  by  8.5  inches  illuminated  area  with  512  electrodes  each  in  both  X  and 
V  directions)  to  generate  two-dimensional  Walsh  function  patterns  and  produce  5 1 2 
by  512  Walsh  coefficients  in  14  seconds.8  The  PISA  produced  successful  results 
for  a  selected  set  of  targets  representing  manmade  cartographic  features  of  the  types 
stated  above. 


*P.I  .  Chen  and  W.W.  SeemuUer,  “Signal  Signatures  of  Topographic  features  Using  Analog  Technology,"  U.S. 
Army  Engineer  Topographic  Laboratories,  Fort  Belvoir,  VA,  LTL-0185,  May  I979.AD-A076  110. 

2p.I-.  Chen  and  W.W.  SeemuUer,  “Application  of  Walsh  Transforms  for  Topographic  feature  Extraction  Using 
a  Sensor  Array  System,"  IEEE  Transactions  on  Instrumentation  and  Measurement,  pp  52-57.  Vol  1M-29. 
March  1980. 

2P.f .  Chen  and  W.W.  Seemuller,  “Detection  of  Signal  Signatures  of  Cartographic  features,”  presented  at  1980 
Army  Science  Conference,  also  published  in  the  Conference  Proceedings,  West  Point.  NY,  June  1980. 

4P.I  .  Chen,  P.W.  Rohde,  and  W.W.  SeemuUer.  “Prototype  Image  Spectrum  Analyzer  (PISA)  for  Cartographic 
feature  l  xtraction,”  U.S.  Army  fnginecr  Topographic  Laboratories,  fort  Belvoir.  VA,  ETL-0204,  October 
1979,  AD-A 080  729. 

5P.I  .  Chen  and  W.W.  Seemuller.  “Signal  Signatures  of  Topographic  f  eatures  Using  Analog  Technology,”  U.S. 
Army  f  ngincer  Topographic  Laboratories,  Fort  Belvoir.  VA,  ETL-0185.May  1979.  AD-A076  110. 

6p.|  .  Chen  and  W.W.  SeemuUer.  “Application  of  Walsh  Transforms  for  Topographic  f  eature  l  xtraction  Using  a 
Sensor  Array  Syster...”  IEEE  Transactions  on  Instrumentation  and  Measurement,  pp  52-57.  Vol  lm-29. 
March  1980. 

7P.l  .  Chen  and  W.W.  Seemuller,  “Detection  of  Signal  Signatures  of  Cartographic  f  eatures.”  presented  at  1980 
Army  Science  Conference,  also  published  in  the  Conference  Proceedings.  West  Point.  NY,  June  1980. 

8P.f  .  Chen.  f'.W.  Rohde,  and  W.W.  .Seemuller.  “Prototype  Image  Spectrum  Analyzer  (PISA)  for  Cartographic 
feature  lxtraction,”  U.S.  Army  Ingincer  Topographic  Laboratories,  lort  Belvoir.  VA.  I  TL-0204,  October 
1979.  AD-A080  729. 


The  sensor  array  minicomputer  system  provides  a  variable  image  threshold  that 
results  in  better  control  of  the  input  images.  The  cartographic  features  as  described 
with  background  scenes  and  noise  were  extracted  successfully  from  aerial  imagery. 
Since  the  sensor  array  minicomputer  system  proved  superior  to  the  PISA  in  this 
application,  a  classifier  was  implemented  for  the  former  to  become  a  semi-automated 
cartographic  feature  extraction  and  recognition  system.9 


In  the  referenced  research,  the  objective  was  to  detect  road  intersections  and 
straight-line  roads  oriented  to  have  0,  45,  90,  and  135  degree  angles  with  respect 
to  the  horizontal  axis  of  the  view  window  only.  In  this  report,  a  refined  scheme  is 
described  to  extend  the  Walsh  transform  processing  technique  to  detect  road  inter¬ 
sections  and  straight-line  roads  at  a  variety  of  angles  with  respect  to  the  axes  of  the 
window.  The  heart  of  the  scheme  was  implemented  as  software  for  an  existing  experi¬ 
mental  system  consisting  of  a  solid  state  sensor  array  as  an  opto-electronic  converter, 
a  minicomputer  as  a  signal  processor,  and  a  computer-controlled  translational  stage 
as  the  imagery  holder.  Successful  extraction  and  recognition  results  are  presented  for 
a  selected  set  of  the  stated  cartographic  features  using  this  experimental  system. 
Finally,  conclusions  are  given. 


SYSTEM  CONFIGURATION 

The  present  system  c  onfiguration  is  shown  in  figure  1.  A  9- by  9-inch  aerial 
transparency,  mounted  on  a  stage  translatable  in  both  the  X  and  Y  directions,  is 
illuminated  by  a  d.c.  -  powered,  white  light  source.  A  section  of  the  transparency  is 
then  imaged  onto  a  Reticon  32-  by  32-element,  solid  state  array  sensor.  The  image  is 
sampled  and  digitized  into  1024  elements  at  1024  gray  levels  and  fed  into  a  Hewlett- 
Packard  (HP)  2108  minicomputer.  A  variable  threshold  is  then  applied  to  the  digitized 
image,  reducing  it  to  a  two-tone  black-and-white  image.  The  white  pixels  are  assigned 
the  value  100,  and  the  black  elements  are  assigned  0.  This  two-valued  image  is  then 
transformed  to  generate  a  two-dimensional  Walsh  transform  coefficient  matrix  com¬ 
posed  of  32  by  32  elements.  The  low  order,  16  by  16  coefficients  are  then  applied 
to  the  software  classification  algorithm,  which  will  be  described  later.  This  algorithm 
classifies  the  image  into  1  of  16  angles  for  straight-line  roads  or  into  1  of  4  angles 
of  road  intersections,  or  “not  recognized”  if  the  image  does  not  fit  either  classification. 


9P.I  .  Chen  and  W.W.  SccmuUcr,  “Detection  of  Signal  Signatures  of  Cartographic  features,”  Presented  at  1980 
Army  Science  Conference,  also  published  in  the  Conference  Proceedings,  West  Point,  NY,  June  1980. 


6 


FIGURE  1.  System  Block  Diagram. 


SOFTWARE  CLASSIFIER 


Two  arrays  arc  passed  to  the  software  classification  subroutine;  the  Walsh  trans¬ 
form  coefficient  matrix  (IXFRM)  and  the  post-threshold  image  matrix.  First,  line 
roads  are  separated  from  road  intersections  to  form  two  independent  large  eases. 


To  differentiate  between  line  roads  and  intersections,  one  must  summarize 
the  first  8  rows  and  8  columns  of  the  transform  matrix,  such  that 


16 

16 

£ 

II 

Of 

i=  l 

3  = 

T  (i,  j) 


where 

T  (i,  j)  is  the  i  th  row  and  j  th  column  element  of  the  Walsh  coefficient  matrix. 


These  summations  are  then  used  as  the  first  step  in  differentiating  between 
line  roads  and  intersections  by  using  the  formula 


C,  -T  ( 1 ,  1) 

R  =  j 

R,  -T(l,  1) 

TO.  l)j 

T(l,  1)| 

If  A  >  0.7,  B  >  0.7.  C’|  >  50,  and  R,  >  50,  then  the  image  is  defined  as  a 
road  intersection. 


This  test  will  recognize  an  intersection  for  all  cases,  except  one.  When  an  inter¬ 
section  is  alined  at  45  degrees  to  the  viewing  axis,  another  test  must  be  made.  From 
the  image  and  Walsh  transform  of  an  intersection  oriented  at  45  degrees,  it  can  be 
observed  that  R(  and  C(  are  both  rather  small  in  magnitude  (both  less  than  50, 
which  was  the  minimum  value  for  an  image  to  be  classified  as  an  intersection  in  the 
last  test  ).  However,  the  sum  of  diagonal  transform  coefficents  [T(l,  1)  +  T  (2,  2)  + 
.  .  .  +  T  (16,  16)]  is  a  reasonably  large  quantity. 


With  the  sum  of  the  first  five  diagonal  coefficients  defined  as 

5 

D  =  £  T  (i,  j). 

i  =  j=  1 
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an  image  is  recognized  as  an  “intersection  at  45  degrees”  if  one  of  the  following 
expressions  is  met: 


If|T(2,3)|>  AND  lT(2’2,?  <  0.3  AND  T  (2, 3)  +  T(l,  1)  >40 

|  T (1,  1)|  | T  ( 1 ,  1) j 

OR 


T  (3,  2)| 

|  T  (2,  3)| 


0.35  AND 


T  (2,  2)  j 
'  T (1, 1)| 


<  0.3 


AND 


T (1 ,  1 )  +  T  (3,  2)  >40 


OR 


D  >  40  AND 


T  (3,  3)  j 
|  T  (1 ,  D| 


>  0.35 


AND 


T  (2,  2) | 
|T(1,  l)j 


<  0.3t 


The  above  tests  in  combination  with  a  minimum  distance  classifier  will  effec¬ 
tively  detect  0  and  45  degree  angles.  However,  another  test  is  necessary,  and  that 
one  is  performed  in  the  image  domain. 

Let 


32 

C  =  2  A  (1 ,  m) 

m  =  1 


32 

D  =  2  A  (32,  n) 

n  =  1 


where  A  (1,  m)  and  A  (32,  n)  with  m,  n  =  1,  2,  ...  ,32  are  the  1st  and  32nd 
rows  of  the  input  pattern  in  the  image  domain.  Increment  m  and  n,  and  compute 
C  and  D  until  C  =  D  =  300.  Then,  if  n  >  m,  the  image  is  at  a  22.55  degree 
angle,  and  if  m  >  n,  it  is  at  67.5  degrees. 


The  algorithm  for  determining  the  angle  of  line  roads  will  be  discussed  next. 
The  image  can  be  divided  into  four  quadrants  (see  figure  2). 
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A(l, 1) 


QUADRANTS 


FIGURE  2.  Input  Image  Array  Divided  Into  4  Quadrants. 


Lei 


16  16 

Mj  =  2  2  A  (i,  j) 

i=l  j=l 

16  32 

M2  =  2  2  A  (i,  j) 

i=  1  j=17 

32  16 

M3  =  2  2  A  (i,  j) 

i  =  1 7  j  =  1 

32  32 

M  =  2  2  A  (i,  j) 

i  =  17  j  =  17 

where  A(i,j)  is  the  input  image  array. 

Several  tests  are  then  applied  to  dv  ^ermine  whether  the  angle  of  the  image  lies 
(a)  between  0  and  90  degrees,  or  (b)  between  90  and  1 80  degrees. 


TEST  1 

If  Mj  =  0,  and  M2,  Mj,  and  M4  do  not  equal  0, 
then  the  angle  is  less  than  90  degrees. 

If  M2  =  0.  and  M,.  Mj,  and  M4  do  not  equal  0, 
then  the  angle  is  between  90  and  1 80  degrees. 

If  M3  =  0,  and  Mt.  M2,  and  M4  do  not  equal  0, 
then  the  angle  is  between  90  and  180  degrees. 

If  M4  =  0.  and  Mj,  M2,  and  Mj  do  not  equal  0, 
then  the  angle  is  less  than  90  degrees. 

If  the  angle  ot  the  image  falls  into  one  of  the  above  categories,  the  routine  then 
proceeds  to  the  Minimum  Distance  Classifier  (MDC).  Otherwise,  it  proceeds  to  TEST 


II 


TEST  2 


(A)  If  Mj  =  M2  =  0,  and  M3  and  M4  do  not  equal  0,  then 

compute  M3  and  M4  until  M3  =  M4  =  1000,  where 

%  is  the  net  increment  of  the  subscript  used  for 
the? outer  summation  to  compute  M3  to  produce 
a  sum  of  1 000. 

*M4  is  similar  to  1^,  but  for  M4 
If  Ij^3  >  lj^4 ,  the  angle  is  less  than  or  equal  to  90  degrees. 
Im3  <  ^M4 » the  angle  is  between  90  and  1 80  degrees. 
Proceed  to  the  MDC  routine. 


(B)  If  M3  =  M4  =  0,  and  M,  and  M2  do  not  equal  0,  then 

compute  Mj  and  M2  until  =  M2  =  1000,  where 
l\f  and  are  similar  to  those  (A)  above. 

If  Im  j  >  *M2>  the  angle  is  less  than  or  equal  to  90  degrees. 

^Mj  ^  the  angle  is  between  90  and  180  degrees. 

Proceed  to  the  MDC  routine. 

If  neither  of  the  conditions  in  (A)  or  (B)  were  true,  then  do  TEST  3. 
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TEST  3 


(A)  If  Mj  =  M3  =  0,  M2  and  M4  do  not  equal  0,  then 

compute  M2  and  M4  until  M2  =  M4  =  1000.  Let 
and  be  the  net  increment  of  the  subscript 
usecl  to  compute  M2  and  M4  respectively. 

If  Jm  >  ,  the  angle  is  less  than  or  equal  to 

90  degrees.  4 

If  <  ^m4>  the  angle  is  between  90  and  180  degrees. 
Proceed  to  MDC  routine. 


(B)  If  M2  =  M4  =  0,  and  Mj  and  M3  do  not  equal  0,  then 

compute  Mj  and  M3  until  Mj  =  M3  =  1000,  where 
Jmi  and  are  similar  to  those  in  (A). 

If  >  Jm3>  the  angle  is  less  than  or  equal  to  90  degrees. 

If  Jjvii  <  Jjvi3  ,  the  angle  is  between  90  and  180  degrees. 

Proceed  to  the  MDC  routine. 


TEST  4 

This  test  is  performed  if  each  of  the  four  quadrants  contains  at  least  one  white 

pixel. 

If  Mj  =  0,  M2  =  0,  M3  =  0,  or  M4  =  0,  then 
proceed  to  TEST  5. 

If  T  (2,  2)  <  0,  the  image  is  less  than  90  degrees. 

If  T  (2,  2)  >  0,  the  image  is  greater  than  90  degrees. 

Proceed  to  the  MDC  routine. 
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TEST  5 


If  Mj  does  not  equal  0,  and  M2  =  M3  =  M4  =  0, 
then  the  angle  is  less  than  90  degrees. 

If  M4  does  not  equal  0,  and  Mj  =  M2  =  M3  =  0, 
then  the  angle  is  less  than  90  degrees. 

If  M2  does  not  equal  0,  and  Mj  =  M3  =  M4  =  0, 
then  the  angle  is  between  90  and  1 80  degrees. 

If  M3  does  not  equal  0,  and  Mj  =  M2  =  M4  =  0, 
then  the  angle  is  between  90  and  1 80  degrees. 

Proceed  to  the  MDC  routine. 

If  the  image  has  not  passed  any  of  the  above  tests,  then  it  is  classified  as  “not 
recognized.” 


THE  MINIMUM  DISTANCE  CLASSIFIERS 
(MDC) 

Separate  minimum  distance  classifiers  (MDC)  were  designed  for  both  linear- 
road  patterns  and  road-intersection  patterns  respectively.  The  set  of  linear-line  roads 
consists  of  16  classes,  and  each  class  represents  a  linear  road  oriented  at  a  particular 
angle  with  respect  to  the  axes  of  the  inspecting  window.  Likewise,  the  set  of  road 
intersections  consisted  of  four  classes.  Assuming  that  these  classes  are  representable 
by  prototype  (or  reference)  patterns  ,  Z2 ,  .  .  .  ,  j£m  ,  then  for  linear  roads, 
m  =  16  and  for  road  intersections,  m  =  4. 

The  Eclidean  distance  between  an  arbitrary  pattern  vectors  x  and  the  i  th 
prototype  is  given  by 
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where 


(  x.  -  )T  is  the  transpose  of  (  x.  -  *  ). 


A  minimum  distance  classifier  computes  the  distance  from  a  pattern  X.  of 
unknown  classification  to  the  prototype  of  each  class  and  assigns  the  pattern  to  the 
class  to  which  it  is  closest.  In  other  words,  X  >s  assigned  to  class  “  i  ”  if  D;  < 
D.  for  all  j  4  i. 


To  establish  the  prototypes  ^ ,  j|2 ,  •  •  • »  ,  one  should  use  the  first  8  rows  and  8 

columns  of  the  summarized  Walsh  transform  coefficients,  Rx,  R2  . . . ,  Rg,  C, ,  C2 
.  .  .  ,  Cg,  as  measurement  vectors.  They  are  expressed  as 

■^i  =  I  ^li»  ^2i>  •  •  '  ’  ^8i>  ^li>  ^2i’  •  •  '  ^8i  1 

Sixteen  prototypes  for  linear  roads  having  incremental  angles  of  1 1.25  degrees  (star¬ 
ting  at  0  degree  with  respect  to  the  horizontal  axis  of  the  window)  were  obtained 
using  the  above  equation.  Likewise,  four  prototypes  for  road  intersections  at  angles 
of  0,  22.5,  45.0,  and  67.5  degrees  with  respect  to  the  horizontal  axis  of  the  win¬ 
dow  were  computed.  These  2.  were  then  used  to  compute  D;  and  classify  the 
incoming  unknown  pattern,  X-  The  complete  software  package  for  this  method  is 
attached  in  the  appendix. 


EXPERIMENTAL  RESULTS 

Targets  representing  linear  roads  and  road  intersections  at  various  angles  with 
respect  to  the  axes  of  the  window  were  constructed  on  a  glass  plate.  The  targets  were 
used  as  input  test  images.  The  glass  plate  was  mounted  on  a  computer-controlled 
translational  stage  so  that  the  relative  position  of  the  images  against  the  surface  of 
the  solid  state  array  can  be  adjusted  arbitrarily. 
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In  figures  3  through  18,  the  original  images  in  spatial  domain,  the  recognition 
results,  and  the  Walsh  transforms  for  linear  roads  at  a  variety  of  orientation  with  the 
window  are  shown.  The  experimental  results  for  road  intersections  having  various 
angles  with  the  axes  of  the  window  are  shown  in  figures  19  through  22.  Correct 
recognition  was  obtained  for  all  cases. 


A  random  translation  of  these  images  against  the  window  was  also  tested.  It 
was  discovered  that  for  the  well-defined  images,  such  as  targets  used  for  this  experi¬ 
mentation,  no  misclassification  was  obtained.  In  figures  23  and  24,  the  typical 
test  results  for  a  linear  road  and  a  road  intersection  translated  away  from  their  center 
position.  Again  correct  recognition  was  obtained  for  both  cases.  Other  tests  for  trans¬ 
lation  were  performed  and  similar  results  were  obtained. 


The  detectable  angular  resolutions  for  the  line  roads  and  the  road  intersections 
were  found  to  be  11.25  and  22.5  degrees,  respectively.  When  the  test  images  were 
corrupted  by  a  large  amount  of  noise,  the  method  failed  to  produce  correct  recognition 
results.  The  amount  of  noise  may  be  reduced  by  thinning  techniques.  However,  the 
thinning  techniques  were  not  tested  during  the  experiments. 


CONCLUSIONS 

1.  Noise-free  linear  cartographic  features,  such  as  roads  and  road  intersections 
at  a  variety  of  angles  with  the  axes  of  the  inspecting  window,  can  be  detected  using 
an  extended  Walsh  transform  processing  technique. 

2.  Minimum  distance  classifiers  together  with  some  pretesting  conditions  were 
found  to  be  very  suitable  for  this  application. 

3.  All  the  test  images  (targets)  at  a  variety  of  angles  and  also  with  arbitrary 
translations  were  detected  and  classified  correctly. 

4.  The  detectable  angular  resolutions  for  linear  roads  and  road  intersections 
were  found  to  be  11.25  and  22.5  degrees,  respectively. 
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FIGURE  3.  Spatial  Signature,  Recognition  Result,  and  Walsh  Transform 

for  Line  Road  at  0  Degree. 


FIGURE  4.  Spatial  Signature,  Recognition  Result,  and  Walsh  Transform 

for  Line  Road  at  1 1 .25  Degrees. 
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FIGURE  6.  Spatial  Signature,  Recognition  Result,  and  Walsh  Transform 

for  Line  Road  at  33.75  Degrees. 
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FIGURE  7.  Spatial  Signature,  Recognition  Result,  and  Walsh  Transform 

for  Line  Road  at  45  Degrees. 
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FIGURE  9.  Spathil  Signature,  Recognition  Result,  and  Walsh  Transform 

for  Line  Road  at  67.5  Degrees. 
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FIGURE  10.  Spatial  Signature,  Recognition  Results,  and  Walsh  Transform 

for  Line  Road  at  78.75  Degrees. 
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FIGURE  1 1.  Spatial  Signature,  Recognition  Result,  and  Walsh  Transform 

for  Line  Road  at  90  Degrees. 
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FIGURE  12.  Spatial  Signature.  Recognition  Result,  and  Walsh  Transform 

for  Line  Road  at  101.25  Degrees. 


WED  JUN  li  198*  14 i 34 


FIGURE  13.  Spatial  Signature,  Recognition  Result,  and  Walsh  Transform 

for  Line  Road  at  112.S  Degrees. 
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FIGURE  14.  Spatial  Signature.  Recognition  Result,  and  Walsh  Transform 

for  Line  Road  at  123.75  Degrees. 
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FIGURE  15.  Spatial  Signature.  Recognition  Results,  and  Walsh  Transform 

for  Line  Road  at  135  Degrees. 
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FIGURE  16.  Spatial  Signature,  Recognition  Result,  and  Walsh  Transform 

for  Line  Road  at  146.25  Degrees. 
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FIGURE  17.  Spatial  Signature,  Recognition  Result,  and  Walsh  Transform 

for  Line  Road  at  157.5  Degrees. 
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FIGURE  19.  Spatial  Signature,  Recognition  Result,  and  Walsh  Transform 

for  Road  Intersection  at  0  Degrees. 


FIGURE  20.  Spatial  Signature,  Recognition  Result,  and  Walsh  Transform 

for  Road  Intersection  at  22.5  Degrees. 
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FIGURE  2! .  Spatial  Signature,  Recognition  Result,  and  Walsh  Transform 

for  Road  Intersection  at  45  Degrees. 
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FIGURE  22.  Spatial  Signature,  Recognition  Result,  and  Walsh  Transform 

for  Road  Intersection  at  67.5  Degrees. 
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FIGURE  23.  Spatial  Signature,  Recognition  Result,  and  Walsh  Transform  for  Line  Road 
at  1 1 .25  Degrees,  and  Translated  to  the  Top  of  Window. 


FIGURE  24.  Spatial  Signature,  Recognition  Result,  and  Walsh  Transform  for  Road  Intersection 

at  0  Degree,  and  Translated  to  the  Right  of  Window. 


APPENDIX 

SUBROUTINE  MATCH 


39 


noon no  non  non 


FTN4*L 


SUBROUTINE  HATCH  < IXFRM.CONS. LUOT * INBUF > 

DIMENSION  IXFRM< 1024 ) * CONS ( 20 >16) » C0LM( 16) > INBUF ( 1024) 

IT0T=0 
DO  20  J=1  *  8 
ISUM=0 

DO  10  1  =  1  »  16 

10  ISUM=ISUM+IABS<IXFRM<  <l-i)*16+J> ) 

C0LH< J)=ISUM 
20  IT0T=IT0T+ISUM 
DO  28  J=9* 16 
ISUM=0 

DO  26  1=1*16 

26  ISUM=ISUM+IABS< IXFRM<  < J-9)*16+I ) > 

C0LM( J)=ISUM 
28  IT0T=IT0T+ISUM 

SEE  IF  INTERSECTION 

A=ABS<  C0LM< 1 )-IXFRM( 1 ) >/IABS( IXFRMC 1 ) ) 

B=ABS( COLM<  9 )-IXFRM< 1 ) )/IABS ( IXFRM( 1 ) ) 

IF  (A.GT.0.7  .AND.  B.GT.0.7  .AND.  COLM( 1 > .GE .50. 

*  .AND.  <C0LM<9) . GE.50) )  GOTO  800 

COMPUTE  DIAGONAL 

IDIAG=0 
DO  900  1=1*5 

900  IDIAG=IDIAG+IXFRM( <I-1)*16+I> 

AT1=IABS(IXFRM(19>) 

AT2=I ABS( IXFRM( 18) ) 

AT3=IABS(IXFRM<34>> 

AT4=IABS( IXFRM<35) ) 

AT5=IABS< IXFRM< 1 ) ) 

IF  (((AT1/AT5) .GT.0.35  .AND.  (AT2/AT5) .LT .0.3  .AND. 

«  ( AT5+AT1 ) .GT.40) 

# » OR .  ((AT3/AT1) .GT.0.35  .AND.  < AT2/AT5) .LT.0.3  .AND.  ( AT5+AT3) 
*. GT.40)  .OR,  UDIAG. GT.40  .AND.  ( AT4/AT5) .GT.0.35  .AND. 
*(AT2/AT5). LT.0.3))  GOTO  810 


COMPUTE  QUADRANTS  TO  DETERMINE  IF  ANGLE  IS  LESS  THAN  90  DEGREES 
OR  IF  IT  IS  BETWEEN  90  AND  180  DEGREES 


M1=0 

M2=0 

DO  120  1=0*15 
DO  100  J=1 *  16 
100  Ml =M1+ INBUF ( I432+J) 
DO  110  J=17 *32 
110  M2=M2+INBUF(I*32+J) 
120  CONTINUE 
M3=0 
M4=0 
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DO  150  1=16*31 
DO  130  J=l*16 
130  M3=M3+INBUF(I*32+J> 
DO  140  J=17»32 
140  M4=M4+INBUF(I*32+J) 
150  CONTINUE 

TEST  ONE 


K=0 


IF 

(M2.NE.0  .AND. 

M3.NE.0  .AND. 

( (Ml. EQ.O 

.AND. 

M4.NE.0) 

.OR. 

«(M4 

.EQ.O  .AND.  Ml 

. NE.O) ) )  K*1 

IF 

(Ml .NE . 0  .AND. 

M4.NE.0  .AND. 

((M2. EQ.O 

.AND. 

M3. NE.O) 

.OR. 

♦(M3.EQ.0  .AND.  H2.NE.0)))  K=8 
IF  <K  .NE.  0)  GOTO  29 

TEST  TWO 

IF  (Ml .NE .0  .OR.  M2.NE.0  .OR.  M3.EQ.0  .OR.  M4.EQ.0)  GOTO  300 

M3*0 

M4=0 

DO  210  11-16*31 
DO  200  Jl*l*16 
M3=M3+INBUF < I 1*32+ J1 ) 

IF  <M3  .GE.  1000)  GOTO  220 
200  CONTINUE 
210  CONTINUE 
220  DO  240  12*16*31 
DO  230  J2=17*32 
M4=M4+INBUF  < I2*32+J2) 

IF  (M4  . GE.  1000)  GOTO  250 
230  CONTINUE 
240  CONTINUE 
250  M3=I 1*32+J1 

M4=I2*32+J2-16 

K=8 

IF  (M3  .GE.  M4 )  K=1 
GOTO  29 

300  IF  (M3.NE.0  .OR.  M4.NE.0  .OR.  M1.EQ.0  .OR.  M2.EQ.0)  GOTO  400 
M1=0 
M2=0 

DO  320  11=0,15 
DO  310  Jl*l *  16 
M1=M1+INBUF ( I 1*32+ J1 ) 

IF  (Ml  .GE.  1000)  GOTO  325 
310  CONTINUE 
320  CONTINUE 
325  DO  340  12*0*15 
DO  330  J2=17 ,32 
M2=M2+INBUF ( I2*32+J2) 

IF  (M2  .GE,  1000)  GOTO  350 
330  CONTINUE 
340  CONTINUE 
350  GOTO  250 
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TEST  THREE 

400  IF  (Ml.NE.O  .OR.  H3.NE.0  .OR.  M2.EG.0  .OR.  M4.EQ.0)  GOTO  500 
M2=0 
M4=*0 

DO  420  11=17.32 
DO  410  J1=0. 15 
M2=M2+INBUF ( J1432+I1 ) 

IF  <M2  .GE.  1000)  GOTO  430 
410  CONTINUE 
420  CONTINUE 
430  DO  450  12=17.32 
DO  440  J2=16»31 
M4=M4+INBUF(J2*32+I2) 

IF  <M4  .GE.  1000)  GOTO  460 
440  CONTINUE 
450  CONTINUE 
460  GOTO  250 

500  IF  (M2. NE.O  .OR.  M4.NE.0  .OR.  M1.EQ.0  .OR.  M3.EQ.0)  GOTO  600 
M1=0 
M3=0 

DO  520  11  =  1 » 16 
DO  510  Jl=0.15 
M1=M1+INBUF(J1*32+I1) 

IF  (Ml  .GE.  1000)  GOTO  530 
510  CONTINUE 
520  CONTINUE 
530  DO  550  12=1.16 
DO  540  J2=16. 31 
M3=M3+ 1 NBUF <  J2*32+ 1 2 ) 

IF  (M3  .GE.  1000)  GOTO  560 
540  CONTINUE 
550  CONTINUE 
560  GOTO  250 

TEST  FOUR 

600  IF  (M1.EQ.0  .OR.  M2.EQ.0  .OR.  M3.EQ.0  .OR.  M4.EQ.0)  GOTO  700 
620  IF  ( IXFRM( 18)  .LE.  0)  K=1 
IF  (IXFRMU8)  .GT.  0)  K=8 
GOTO  29 


TEST  FIVE 


700  IF 

(M2.EQ.0  .AND.  M3.EQ.0 

.AND. 

( (Ml.NE.O 

.AND. 

M4.EQ.0) 

.OR 

*(M4 

. NE.O  .AND.  Ml.EG.O))) 

K=1 

IF 

(M1.EQ.0  .AND.  M4.EQ.0 

.AND. 

( (M2. NE.O 

.AND. 

M3.EQ.0) 

.OR 

*(M3 

.NE.O  .AND.  M2.EQ.0) ) ) 

K=8 

IF  (K  .NE.  0)  GOTO  29 
WRITE  (1.710) 

710  FORMAT  C  IMAGE  DOES  NOT  FIT  ANY  TESTS') 
RETURN 

ROAD  INTERSECTIONS 


non 


810  MINNUM=19 
GOTO  65 
800  K=17 

FIND  CLOSEST  MATCH 

29  AMIN=1000. 

L=3 

IF  (K.EQ.l  .OR.  K.EQ.8)  L*8 

DO  40  I=K*K+L 

DIF*0. 

IF  (I  .EQ.  19)  GOTO  40 
DO  30  J=1 > 16 

30  DIF=DIF+(C0LM(J)/IT0T-C0NS(I»J))**2 
IF  (DIF  .GE.  AMIN)  GOTO  40 
AMIN=DIF 

MINNUM=I 
40  CONTINUE 

IF  (MINNUM  .GT.  16)  GOTO  60 

IF  (MINNUM  .LT.  8)  ANGLE=MINNUM*11 .25 

IF  (MINNUM  .EQ.  8)  ANGLE-0. 

IF  (MINNUM  .GT.  8)  ANGLE=(MINNUM-1)*11.25 

WRITE  (LU0T.50)  ANGLE 

50  FORMAT  C  LINE  ROAD  AT  *»F6.2»*  DEGREES*) 
RETURN 

60  IF  (MINNUM. EQ. 17)  GOTO  65 
C=0 
D=0 

DO  62  1=1 i 32 
C=C+INBUF ( I ) 

D=D+INBUF (992+1) 

IF  (C.EQ.300  .OR.  D.EQ.300)  GOTO  63 

62  CONTINUE 

63  MINNUM=18 

IF  (D  .GT,  C)  MINNUM=20 
65  ANGLE* ( MINNUM- 1 7  >  #22 . 5 
WRITE  (LUOT t70)  ANGLE 

70  FORMAT  C  INTERSECTION  AT  *.F6.2»*  DEGREES*) 
RETURN 
END 
END$ 


43 


